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Abstract 

The relationship between the parton transverse momentum and 
'!::j- ■ longitudinal momentum, which is obtained by analysing the hadron 

^ . tensor in the deep inelastic scattering (DIS) process, indicates that the 

transverse size of partons in a nucleon is sensitively affected by the ef- 
0> ' fective nucleon mass. The change in this transverse size gives rise to an 

r~| ' Additional Parton Evolution (APE) in the nuclear environment refer 

Qh! to the normal QCD parton evolution thus resulting in new parton dis- 

Q^' tributions. The self-consistent integral equations for these new parton 

^ ■ distributions are given. Having undergone the normal QCD evolution 

and the additional evolution in the nuclear environment, the partons 
with the small x participate in recombination and contribute to the 
^ ■ nuclear shadowing effect. The reasonably good descriptions of the nu- 

clear effect of structure functions, nuclear Drell-Yan ratios and nuclear 
gluon distributions are presented. A quantitative comparison between 
the APE model and the Q^-rescaling model is made. 
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1 Introduction 

Since the discovery of the nuclear modification in the nucleon structure func- 
tion F2{x, Q"^) by the European Muon Collaboration (EMC) m, various mod- 
els have been proposed to investigate the nuclear effect. Several years ago, the 
New Muon Collaboration (NMC) accurately measured the structure func- 
tions F2 at very small x. These data as well as the EMC data 0] provide an 
insight for investigating the physics details in the shadowing region(x < 0.1). 
The need for testing the models of the EMC effect becomes more and more 
important for further understanding of the nature of the nuclear effect. Bick- 
erstaff et al. |^ found that although most of the theoretical models can 
provide good explanations for the EMC effect, they can not describe the 
nuclear Drell-Yan process well. The NMC [H analyzed the inelastic J/ip 
production data and obtained the gluon distribution ratios Gsn{x)/Gc{x) 
which indicate a modification of the gluon distribution in the nuclear envi- 
ronment. Kumano has made progress towards the unified description of 
the EMC effect and the nuclear gluon distribution in the Q^-rescaling model 
with the recombination effect. Although the Q^-rescaling model is a useful 
and effective model in explaining deep-inelastic data in the valence domi- 
nated medium x region, it is not clear whether the momentum distributions 
of sea quarks and gluons follow the same rescaling as that for valence quarks 
in the nuclear environment. According to the QCD theory, sea quarks have 
different dynamic properties while interacting with gluons. There is a con- 
straint condition on the number conservation for valence quarks but not for 
sea quarks. So it is reasonable to believe that the modification of sea quark 
distributions should be different from that of valence quarks. For further 
discussion, we assume that the transverse size of partons in a nucleon can 
be described by an average parton transverse momentum which is related to 
the longitudinal distributions of partons. The relationship between the par- 
ton transverse momentum and longitudinal momentum, which is obtained 
by analysing the hadron tensor in the deep inelastic scattering process, indi- 
cates that the transverse size of partons in a nucleon depends sensitively on 
the effective nucleon mass. The change of this transverse size gives rise to 
an additional parton evolution in the nuclear environment refer to the nor- 
mal QCD parton evolution thus resulting in new parton distributions. The 
probabilities of the additional evolution are directly related to the ratio of 
the average parton transverse momentum in the bound nucleon to that in 



the free nucleoli. According to these probabihties, we give the self-consistent 
integral equations for parton distribution modifications due to the additional 
evolution in the nuclear environment. 

Considering the fact that the nucleon diameter is about 1.8/m, in an 
infinite momentum frame, one can estimate that the longitudinal size of a 
nucleon in a Lorentz contracted nucleus is 

D = 1.8{fm)mN/pN, 

and the longitudinal localization size of a parton with the momentum xpN 
is L = 1/{xpn)- If a parton has a small x and consequently its dimension 
(L) exceeds the nucleon longitudinal size (D), it leaks out of the nucleon and 
recombines with partons of other nucleons. We suppose that the partons 
participating in recombination have undergone the normal QCD evolution 
and the additional evolution in the nuclear environment before the fusion. 
After considering the additional evolution in the nuclear environment with 
the parton recombination effects, we find that the new model can give a well 
explaination to the nuclear effects of the structure function, nuclear Drell- 
Yan process and inelastic J/ip production. Finally, a quantitative comparison 
between the APE model and the Q^-rescaling model is made. 

2 Additional Parton Evolution(APE) Due to 
Nuclear Environment 

Partons in a free nucleon have their own distributions with a specified scale. 
It is generally accepted that the transverse size of partons in a nucleon is 
affected by the nuclear environment. We introduce an average transverse 
momentum to describe the change of the transverse size of partons. 

In order to obtain the average transverse momentum of partons, let us 
repeat the ^ variable analysis |^ of the hadron tensor M/^,^ in the deep inelastic 
scattering process shown in Fig.l . W^i, can be written as 



W,, = ^Jd'xe"^-^p\[J^{x),3M]\p) 



-9,. + ^)W^i + -V(P. - ^^.)(P. - ^I^W. (1) 

q iiij^ q q 



where rriN is the mass of a nucleoli. If we take the frame in which p = {rriN, 0) 
and neglect the quark mass, then 

W,.^ f^fi^)6iq' + 2k-q)w,. (2) 

where 

w^u = [{k + q)^,K + {k + q)yk^ + —g^,y\ (3) 

is the quark tensor in the lab. frame and f{k ■ p) is a Lorentz-invariant wave 
function. To extract Wi and W2 from Wfj,^, we introduce two contractions 

B = g'"'W^, = -3Wi + (l + ^)W2 (4) 

and then Wi and W2 can be expressed as 

W,= ^-iC^B) (6) 

W2 = ^^"^ (7) 

' 2(1 + z/Vg2) V ) 

If we take the frame in which q = (z/, 0, 0, ^/^u^~+~Q^) and A; = (/cq, 0, ko sin 6', /cq cos i 
then 

By means of Eq.(|^) and the contractions defined by Eq.(^ and Eq.(|^), one 
can find the relation 

(*i) - -^ (9) 



Taking Wi and W2 in Eq.(|) and Eq.(||) as those in Refs.0 and |T0[ and 
neglecting 0{m%/Q^) and OiQ"^ /v"^) terms, one easily obtains 



Furthermore, we assume that there is a similar relationship between the gluon 
transverse momentum and longitudinal momentum distributions, i.e.. 



nN 



Eq.(|TD|) and Eq.([n]) with the bound nucleon effective mass mj^(^A) ("^Af(i) = 
771 at) have been written in a unified form for partons in both a free nucleon(N) 
and a bound nucleon(A) in a nucleus with the mass number A. The effective 
mass of the bound nucleon is less than that of the free one, i.e., in terms 
of Eq. (0) and Eq. (|ll|) , the average transverse momentum of partons in the 
bound nucleon is less than that in the free nucleon, which mainly causes the 
parton evolution in the nuclear environment. By using the same technique 
in deriving the famous Altarelli-Parisi evolution [^, we obtain the evolution 
probabilities for the evolution from a parton with momentum fraction y to 
another parton with the momentum fraction x; 

^''-'^^y^ - ^^'^""^^a^W ~'^ " ~y^^ {Kiix,Q^)rt ^''^ 
fg^a-) = —[{-? + (1 - ^)2]/^Mi^i^!))!^ (13) 

X a 1 + (1- ^)2 (kVx 0^)f 

Fg^gC-) = -C^2(G)[-i- + ^i^ + -(1 - -) 

y T^ (^~y)+ y y y 

where the strong interaction coupling constant a^ is 

a,(g2)=47r/[/?o/n(gVA2)] 

with A = 0.19Ge\/ and /5o = 11 - 2Nf/3,Nf = 3. C2{R) = | and C2{G) = 3 
are color factors for A'"c = 3. Therefore, the quark distributions modified 



by the additional parton evolution(APE) in the nuclear environment can be 
written as follows: 



qr^ix,Q') = q:{x,Q') + C F,^^,S-fl^^^^^dy (16) 

Jx y y 

qr^{x,Q^) = q^{x,Q^) + /'f,^,(-)^2^^^^ 

Jx y y 

+ i Fo-.,^(-)^-^dy (17) 

, Jx V V 



fi X G%j^ 
+ / Fg^g{-) dy (18) 

J^ y y 

where the input parton distributions q^{x, Q"^) and G^{x, Q"^) are taken from 
Refs. [|12[- [|l^]. In fact, these equations are the self-consistent integral equa- 



tions since the evolution probabilities Fp.^p^ (z) contain the parton transverse 
momenta which are related to the parton longitudinal momentum distribu- 
tions g^(^)(x,g2) and G^^^\x,Q^) by Eq.(|lOD and Eq.(|llD. In the calcula- 
tion of the parton distributions in the bound nucleon, the numerical solution 
of these equations is obtained by the iteration method. 

3 Leak-out Sea Quarks and Gluons in Nu- 
clear Environment 

Although the leakage in the nuclear environment can occur for all partons, 
the most important contribution arises from partons with the largest spatial 
uncertainty, i.e., those with the small x. We assume that the leak-out partons 
are sea quarks and gluons for which the momentum cutoff function |jl5| is 
taken as 

(3{x) = exp(-^^i|i^), 
namely, the distributions of leak-out (LK) partons are 

/^(a;,Q^)=/5(x)/^^(x,Qg) 



where the input parton distributions p^^^ = uf^^, u"^^^, df^^, d"^^^, s^^^ 
, s^^^ and G^^^ together with u^^^ and d^^^ are those having undergone 
the normal QCD parton evolution and the additional parton evolution due 
to the nuclear environment. The momentum cutoff function is shown as a 
function of x in Fig.2 with various Hq. As discussed in section 1 , the partons 
with X < 0.11 can leak out of a nucleon. We then find in Fig.2 that an 
appropriate choice of j/o is 2.0/m. 

4 Parton Recombination and Further Mod- 
ification to Parton Distributions 

According to the original idea of the parton recombination by Close, Qiu 



and Roberts ||T5|, a leak-out parton can fuse with a parton from another 
nucleon. In general, the modification of a parton distribution p^{x^), due to 
the process of producing the parton p^ with the momentum X3 by fusion of 
partons pi and ^2, is given by [^ 



^P?.{x-i) = K J dxidx2Pi{xi)p2{x2) 



xTp^p2^p^{Xi,X2,X3 = Xi +X2)S(X3,XI,X2), 

where K is given by Refs. |jl5| and ||T6i , 



K = 9A'/'as/i2Ro'Ql) 
The nuclear radius is i? = RqA^'^ with Rq = 1.1 fm |l^. The function 



S{xs,xi,X2) defined by Ref. p is introduced to keep the momentum con- 
servation. The parton fusion function Tp-^^p^^p-^^xi, X2, X3) is a probability for 
producing a parton ^3 with momentum X3 by a fusion of partons pi and p2 
with momentum Xi and X2, respectively. It is related to a splitting function 
Ppj^^p^{z) in the Altarelli-Parisi equations |jTll by 

where Cp^p^^p^ is the ratio of color factors in the processe piP2 -^ Pz- Using 
the fusion functions, one can obtain the modifications of parton distributions 



due to the recombination(RC)0. With our notations, the exphcit expres- 
sions for the modifications /S.q^^{x) and AG^^{x) are 

6 JO x' 

_ I 2 

+G^^^(x')g^'^(x-x')][l+ " 



X 



^ /'^^a;a;'[/^i^(a;)G'^i^(^') + g^^(a;)G'^^^(a;')]^_[i + [. ^ I'l 



6 JO x' X + x' X + x' 

^^x /'dx'xx'[g^^^(x)g^^(x') + g'^'^(x)g^^^(x')]^^^^ (20) 



9 JO (x + x') 

and 

xAG'^^(x) 

= ^xj dxV(x-x')[G^^''(x')G^'^(x-x') 

,,, 1 |. x' X — x' x'(x — x'). 









^^x /'rfx'xx'[G^^(x)G'^^^(x') + G^''''(x)G'^^(x')] 



4 Jo 

Xl ^ — + — + 



(x + x')^ x' X (x + x')^ 

+ — X / rfx'y[xV''''(a;')(a;-a;')g'''^(a;-x') 
9 7d Y 

+x'qf''''{x'){x - x')qf''ix - x')]^[x" + (x - x')'] 



K fi 
"6 

K /•! 
TJo 



/' rfx'xG^^^(x) y^ix'qt^'ix') + x'gf^(x')]^-[l + [^^] 
dx'xG^^(x) ^[x'g,^^^(x') + x'g^^^(x')] 



1 ' 2 

1 r. , X 



X- 



_[l + r^^ii (21) 

x + x'^ ^x + x'^ ^ ^ 

The exphcit Q^ dependence in the parton distributions is not shown in the 
above two equations in order to simphfy the notation. Therefore, we write. 



for the nuclear gluon distribution, 

xG^{x, Ql) = xG^^^(x, Ql) + AxG^^(x, Ql) (22) 

Similarly, for the nuclear quark distribution, 

xg^(x, Ql) = xg^^^(x, Ql) + Axg^^(x, Ql) + Axgf ^(x, Ql) (23) 

where Axgf "^(x, Qo) is the modification due to the gluon shadowing(GS) and 
is derived from the evolution equations. A crude estimate of Axgf '^(x, Ql) 



is given by 15 



x*(x.QS) = -^^!^^^ (24) 

Using this relation, we expect that the sea-quark distribution is also affected 
by the gluon modification as follows: 

xAgf ^(x, Ql) = -9{xo - ^^^ ^l^^^^i^^Ql)] ^^5) 

where the step function 9{xq — x) is introduced because the relation is valid 
only at very small x (x less than xq with xo=0.1). 



5 Explanation of Nuclear Effect in Nucleon 
Parton Distributions 

We have considered the additional parton evolution and parton recombina- 
tion in sections 2 and 4, respectively, and gotten the nuclear parton distribu- 
tions at Ql = A.OGeV^. To compare the results with the experimental data 
or display the Q^-dependence of the model, these distributions are evolved 
to those at lager Q^ by using the ordinary Altarelli-Parisi equations [|11[]. For 



the valence quark(g^), sea quark(gs), and gluon(G), the evolution equations 



are 



9 Af .^ [^ dy A, .Nn r^i 



9i<i^(^^t) = l -f^^iy^t)P^^l^] (26) 

|^s^(^,t) = £ %^{y,t)P,,[-] + G^{y,t)P,G[-]] (27) 



^G\x,t) = f %f{y.t)PGA + G\y,t)Pac[l]] (28) 

ot Jx y y y 

where Pij is the sphtting function, and t is defined by 

t = -{2/(3o)H(^s{Q')/as{Ql)] (29) 

We used the solution of these integral equations in the leading order and the 



evolution subroutine of Ref. |T8ll in our numerical calculation. 



5.1 Nuclear Structure Function ^2(0;) 

We now proceed to compute the modification in the structure function. In 
order to compare with the experimental data , let us define the ratio of the 
average nuclear structure function to the deuteron structure function as [|19 

with the average nuclear structure function: 

F2^ix,Q') = j[F,^{x,Q') - ^iN-Z)iF^{x,Q')-Fi{x,Q'))] (31) 

where F^ and F2 are the free neutron and proton structure functions, re- 
spectively. F2a{x,Q'^) is the nuclear structure function. The second term 
compensates for the neutron excess. By considering the Fermi motion of 
nucleons in the nucleus, F2a{x, Q"^) can be written as [T^ 

i^2A(x,g^) = E/^I^A(p-)^Ff^^(f,g^), (32) 

with z = (po + P3)/^TT'N,Po = f^N + ^A , ^A is the separation energy of a 
nucleon in the single-particle state A, and ipxip) is the single-particle wave 
function of the nucleon in the momentum space, which satisfies the light-cone 
normalization: 

"^'^ \MP)fz = l (33) 



(27r)3 
In the following calculation, eA and ipxip) are taken from Ref. |20|- By 



using the nuclear parton distributions after considering the additional parton 
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evolution, recombination and the ordinary Altarelli-Parisi evolution to higher 
Q^, the bound nucleon structure function F2 (x, Q^) can be expressed as 

+(i^(x, Q2)] + 2[s'^(x, Q2) + 5^(a;, q2)]| (34) 

We now calculate R^/^{x^Q'^) in Eg. (pO]) . The deuteron structure function 
F^lXjQ"^) is taken from Ref.||21||. Similar to Ref . ||22|| , we take 



^^(^) = l + 0.481n(2-A-V3] (35) 



where the constant 0.48 is determined by 171^(^00) = 0.75mN p3 . 

The theoretical result in the range 0.1 < x < 0.2 can not be well de- 
scribed since the gluonic modification due to Eq.( 541) is valid at very small 



X. However, using the KMRS-B0[112|, our calculation results indicate that 



the modified parton distribution can be well parametrized by the analytical 
form of the corresponding input distribution. The parameters in the analyt- 
ical form of the modified parton distributions can be obtained by fitting the 
numerical results. So we obtain a smooth curve of the ratio -R*^"/^(x, Q^) 
which is compared with the experimental data of the EMC effect |I[] in Fig.3. 
The result is in qualitative agreement with the experimental data. To show 
the distinct effect of the additional parton evolution, the following results are 
no longer presented by fitting numerical results and taking into account the 
Fermi-motion correction. 

In order to investigate the Q^-dependence of the ratio R^°'^^{x, Q^), the 
results obtained by using the GRV inputs[|l3l and the KMRS-BO inputs [|T2 



are plotted in Fig. 4(a) and Fig.4(b), respectively. One can easily find that 
the result obtained by using the KMRS-BO inputs is less dependent on Q^. 
In Fig. 5, the results with various input parton distributions of different pa- 
rameterizations are presented, and they show that, in the small x region, 
the recombination results after undergoing the additional parton evolution 
are very sensitive to the input sea-quark and gluon distributions. Using the 
Duke-0wens(2) inputs of Ref. iH] with hard gluon distributions, theoretical 
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results underestimate the shadowing, which is not difficult to be understood 



since the parton parameterizations of Ref. |14] did not ffi to the experimental 



data in the small x region. This is the reason why we take KMRS-BO []T2 



or GRV [|T3| parameterizations as input parton distributions in most of the 



following calculations. 



5.2 Nuclear Drell-Yan Ratio 



Several years ago, the E772 Collaboration [Q at Fermilab published the 
data of high-mass dilepton production measured in the nuclear Drell-Yan 
(DY) process. These data aroused special attention in clarifying the different 
explanations for the nuclear effect on the parton distributions. Their results 
show that the ratio of DY dimuon yield per nucleon on a nuclear target to that 
on a free nucleon is slightly less than unity if the momentum fraction x carried 
by a target quark is less than 0.1. The ratios over the range 0.1 < x < 0.3, 
however, do not reveal distinct nuclear dependence. 

In order to show the prediction of nuclear Drell-Yan ratio in our present 
model, let us begin with the definition of the nuclear Drell-Yan ratio : 

rpA/Df \ _ j dx a a [x , x)/dx ax , , 

^^^~ Jdx'd^a^-D{x',x)/dx'dx' ^ ' 

where the differential cross section of the nuclear Drell-Yan process is: 

Hh-A{x,x) (37) 



dx'dx 3 Sx'xM'^ 

with 

H,.Aix', x) = Y: e^x'gf (x')xgf (x) + x'q^ix')xqf{x)] (38) 

i 

In a more specified case for the proton-nucleus reaction: 

Hp_Aix',x,Q^) = j-x'[Au^ix',Q^)+d^ix',Q')]xS''ix,Q^) 
54 

+ ^x'5^(a;', Q')x[{A + 3Z)u^{x, Q') 
+ {AA - 3Z)d^{x, Q2) + 2AS^ix, g')] (39) 



12 



where the total sea quark is S = Ug + Us + ds + ds + s + s. By using the nuclear 
parton distributions in our model, we calculated the nuclear Drell-Yan ratios 
for '^^Ca. The integral range for x' in Eq.(PBD is determined according to 
the kinematic region of the experiment in Ref.p^, i.e. x' — x > 0, and 



0.025 < X < 0.30. To avoid the uncertainty derived by the crude estimate in 
Eq. ( ^4D for taking into account the gluonic modification, the results, except 
those in the range 0.1 < a; < 0.2, are compared with the E772 [^ data in 
Fig. 6. Fig.6 indicates that the APE model as well as the Q^-rescaling model 
can explain the Drell-Yan ratio reasonably well. However, in the small x 
region(x < 0.1), the result obtained by the APE model seems to be better 
than that of the Q^-rescaling model. 

5.3 Nuclear Gluon Distribution 

In 1992, the NMC [Q analyzed the inelastic J/ip production result with 
RiniSn/C) = 1.13 ± 0.08 by the color singlet(CS) model @ and obtained 



the gluon distribution ratios Gsn{.x)/Gc{x). The result sheds light on gluon 
distributions in nuclei. 

The results of Gsn{.x)/Gc{x) calculated in our model together with the 
experimental data from NMC are presented in Fig. 7, which are qualitatively 
similar to those given by the Q^-rescaling model in Ref. Q]. Due to the 
gluon shadowing, the ratios Gsn{,x)/Gc{x) are less than 1 at x < 0.09. It is 
noteworthy that experimental errors are very large in comparison with typical 
theoretical modifications. In order to test our model, we should wait for 
more accurate measurements of Ga{x), for example, a proposed experiment 



at RHIC ||2^ in the small x region for investigating details of the gluon 
shadowing. 

6 Discussion and Summary 

It is worthwhile to compare the results of the present model with those of 
the Q^-rescaling model. Although these two models seem to have the same 
physics origin, they have different descriptions of the nuclear effect on sea 
quarks and gluons. In the APE model, the evolution of valence quarks, 
sea quarks and gluons are determined by different interaction vertex factors, 
and their effects on their own momentum distributions are also different. 
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However, the Q^-rescaling model, ignoring the difference in additional parton 
evolution modes for different kinds of partons, describes the influence of the 
nuclear environment on the distribution functions of valence and sea quarks 
by using the same Q^-rescaling mechanism. In the medium x region where 
the valence quarks dominate, the Q^-rescaling model can explain the EMC 
effect reasonably well, indicating that the nuclear effect on the valence quark 
momentum distribution can indeed be well described by the Q^-rescaling 
mechanism. In the small x region where sea quarks and gluons dominate, 
it is not clear whether the Q^-rescaling picture could be used also for sea 
quarks and gluons. To avoid this uncertainty, we attempt to find a unified 
description of the nuclear effect on sea quarks and gluons as well as valence 
quarks. This is the purpose of our investigation. To make a quantitative 
comparison between the APE model and the Q^-rescaling model, some results 
of these two models are shown in Figs. 8-9 in different x scales. Firstly, Figs. 8- 
9 indicate that both models give reasonable explanations of the EMC effect 
in the medium x region. However, the expectations of the ratio R^"-!^ for two 
different input distributions (GRV and Duke-Owens) are almost the same in 
the APE model but distinctly different in the Q^-rescaling model in which 
the ratio is slightly sensitive to the choice of free parton distributions. In the 
lager x region, there are distinct modifications of the recombination effect in 
the ratio of R^'^i^ which are similar to, but different from the contributions 
of the nuclear Fermi motion (cf. Fig.3). For the small x region, we are 
always seeking for some mechanisms which account for the nuclear shadowing 
effect. Unfortunately, neither the APE model nor the Q^-rescaling model 
can give suitable explanation of the nuclear shadowing effect. It has been 
emphasized by Ref.^ that the shadowing due to the recombinations in the 
nuclear environment is mainly produced as a result of the modification in 
the gluon distribution. Considering this recombination effect, both the Q^- 
rescaling model (see Ref. 0) and the APE model (see Fig.4) can fairly 
explain the EMC and NMC shadowing data by using the KMRS-BO input 
distributions. But, by using harder input gluon distributions such as those 
of the Duke-0wens(2)|T^, the theoretical results (see Figs. 8-9) of these two 
models underestimate the shadowing. It seems, however, that by considering 
the recombination effect, the prediction of nuclear shadowing in the APE 
model is closer to the experimental data than that in the Q^-rescaling model 
with these harder input gluon distributions. This is possibly attributed to 
the better descriptions of the nuclear effect on sea quarks and gluons in our 
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model. 

As for the prediction of the nuclear Drell-Yan ratio, the result in Fig. 6 
indicates that, with the same input parton distributions of GRV, the APE 
model has almost the same prediction of the nuclear Drell-Yan ratio for "^^Ca 
as the Q^-rescaling model in the range 0.2 < x < 0.35 although there is a 
slight difference in the nuclear shadowing region {x < 0.1). The results of the 
nuclear Drell-Yan ratios for other nuclei such as C, Fe, W are not presented 
since the results are qualitatively similar. The results of these nuclear Drell- 
Yan ratios do not show a distinct nuclear dependence, which is consistent 
with the experimental fact. 

In summary, we have put forward the APE model which attempts to 
explain the nuclear effect on the parton distributions by the change of the 
transverse momentum of partons in the nuclear environment instead of the 
Q^-rescaling model. By means of the APE model with the parton recombina- 
tion effect, we investigated the nuclear effect including the nuclear structure 
function F.^{x, Q^) in the whole x region, nuclear Drell-Yan process and nu- 
clear gluon distribution. The results of the nuclear structure function , in the 
small X region {x < 0.06), are of a little Q^-dependence but are very sensitive 
to the input sea-quark and gluon distributions. By using the proper input 
parton distributions, we obtained reasonably good agreements with the ex- 
perimental data of the EMC effect, nuclear Drell-Yan ratio and nuclear gluon 
distributions. 

In addition, the comparison between the APE model and the Q^-rescaling 
model shows, on one hand, that the results are qualitatively similar since they 
have the same physics origin; and on the other hand, with better description 
of the nuclear effect on sea quarks and gluons, the value of R^''^ in the APE 
model is less sensitive to the choice of the free parton distributions than that 
of the Q^-rescaling model. Our present investigation should be considered as 
the first step for seeking for more subtle mechanism with the same physics 
origin as the Q^-rescaling model but with less uncertainties. Further, we are 
looking forward to a unified description of the nuclear effect and the nuclear 
shadowing within the framework of the APE model. 
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Figure Captions 

Fig.l. Deep inelastic scattering in the free parton model. 



9 9 9 



Fig. 2. Momentum cutoff for leak-out partons, (3{x) = exp{- 

Fig.3. Comparisons with SLAG data for '^°Ca with Ferim motion correction. 



Input parton distributions are those of KMRS-BO |T^. Q^ = A.OGeV 
Fig.4(a-b). Comparisons of the ratio R^"/^ with EMC-90 § and NMC 



data by using (a) the GRV [T^ input and (b) the KMRS-BO [|T2| input. 



Q'^ = A.OGeV"^ and S.OGeV^ for the solid and dashed curves, respectively. 

Fig.5. Comparisons of the ratio R^"/^ with EMC-90 § and NMC § data by 
using different input distributions. The solid, dashed and dotted curves are 



for the KMRS-BO Q ,the GRV [Tl and Duke-0wens(2) |T|, respectively. 
Q2 = A.OGeV^. 

Fig. 6. The nuclear Drell-Yan ratios T^'^^^{x) predicted by the present model(solid 
curve) and the Q^-rescaling model (dashed curve) are plotted versus x by us- 
ing the GRV input distributions at Q^=20 GeV^. The experimental data are 
taken from the E772 Collaboration |24 



Fig. 7. Comparisons of the result of Gsn{x)/Gc{x) with NMC data [^ by 
using the KMRS-BO input distributions at Q^=4: GeV^. 

Fig.8(a-b). The ratio R^^-Z^^x) given by the APE model is plotted, (a) in 
the linear x scale; (b) in the logarithmic x scale. Results by using the GRV 
inputs(solid curve) and the Duke-0wens(2) inputs(dashed curve) at Q^=4 
GeV^. 

Fig.9(a-b). The ratio R'-'°'^^(x) given by the Q^-rescaling model is plotted, 
(a) in the linear x scale; (b) in the logarithmic x scale. Results by using 
the GRV inputs(sohd curve) and the Duke-0wens(2) inputs(dashed curve) 
at Q2=4 GeV\ 
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